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Abstract 

The equations describing planar magnetoacoustic waves of permanent form in a cold plasma are 
rewritten so as to highlight the presence of a naturally small parameter equal to the ratio of the 
electron and ion masses. If the magnetic field is not nearly perpendicular to the direction of wave 
propagation, this allows us to use a multiple-scale expansion to demonstrate the existence and nature 
of nonlinear wave solutions. Such solutions are found to have a rapid oscillation of constant amplitude 
superimposed on the underlying large-scale variation. The approximate equations for the large-scale 
variation are obtained by making an adiabatic approximation and in one limit, new explicit solitary 
pulse solutions are found. In the case of a perpendicular magnetic field, conditions for the existence 
of solitary pulses are derived. Our results are consistent with earlier studies which were restricted to 
waves having a velocity close to that of long- wavelength linear magnetoacoustic waves. 

1 Introduction 

For a plasma composed of cold electrons and a single species of cold ions, both collisions and Lan- 
dau damping ca.n be neglected with the result that a two-fluid model provides an accurate description 
i Kakutani et alJll968l) . Such a model is governed by the continuity and momentum equations for elec- 



trons and ions, and Maxwell's equations. In the study of non-relativistic hydromagnetic waves with a 
frequency much less than the plasma frequency, these equations may be simplified somewhat by neglect- 
ing the displacement c urrent and taking the number densities of electrons and ions to be equal, except 
in Poisson's equation l|Kakutani et al.lll967j) . Then taking all quantitie s to be indepen dent of y and z 
one arrives at a set of equations governing planar hydromagnetic waves. Ill'ichevI l)l996|) integrates these 
to obtain the following equations for a magnetoacoustic wave of permanent form propagating in the 
x-direction at a constant speed V: 

Av Ri cose 

nw — KiBz, (la) 



de V 

dw Ri cos 9 
" V 

dBy „ Rp cos 9 



RiBy + R, sin 9 (l-n), (lb) 



R^nw^ — nB^, (Ic) 

d^ V 

— -ReTiv — nBy, (Id) 



where ^ — x — Vt, 



1 1 



„ l-^[B'y+2BySin9 + B!), (2) 

By=By- sin 9, (3) 

and 9 is the angle between the equilibrium magnetic field and the x-axis. In the above, ^, the ion density 
n, the y and z components of the ion drift velocity v and w, and the magnetic field {B^, By, B^) are 
normalized, respectively, by the characteristic length /, the equilibrium ion density, the Alfven velocity 
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Va, and the equilibrium magnetic field strength. In such units the speed of a linear long- wavelength 
magnetoacoustic wave is unity. The remaining parameters are defined by Ri — ujcil/VA and Re — lUccI/Va 
where LUd and Uce are the ion and electron cyclotron frequencies, respectively. The values of the dependent 
variables in the absence of a wave are zero for v, w, By, and B^, and unity in the case of n. 

The results we present here originate from our observation that equations Hla|l may be rewritten so 
as to include a small parameter, e. After introducing the variable s = Re£,/V, where V = VsecO, the 
equations take the form 



dt; 

d^ 
dw; 

d7 

dBy 

ds 
dB, 

ds 



-£{nw + VB^}, (4a) 

e{nv + VBy + V tan 6* (1 - n)}, (4b) 

n{Vw + B^}, (4c) 

-n{Vv + By}, (4d) 



where e = Ri/Re is simply the ratio of the electron and ion masses. With e ~ lO^'^ or smaller, any 
analytical treatment should make use of the smallness of this parameter and some form of perturbation 
theory is obviously called for. Furthermore, even a purely numerical method of solution should attempt 
to take advantage of the smallness of e as otherwise it means integrating over many small-scale variations 
before the underlying variation on the large scale takes place. 

From H4a() it can be seen that generally v and w chang e on a m uch slower scale than By and B^ . It is 
therefore natural to make an adiabatic approximation fsee lHakenI l(l983i) ) which is equivalent to equating 
the right-hand sides of l|lc|) and H4d|l to zero. This gives 

B, = -Vw, By = -Vv. (5) 

These relationships can be used to eliminate By and Bz and it is shown in Appendix that all the 
dependent variables can be expressed in terms of a single variable which itself satisfies a Newtonian type 
energy equation with a polynomial Sagdeev potential. This equation is then used to show the existence 
of solitary pulses and nonlinear waves on the larger scale. However, it should be stressed that the very 
nature of the adiabatic approximation as used in the above is to eliminate any variation on the smaller 
scale. Furthermore, the class of solutions to (|4a|l is restricted, in view of relationships (0), which implies 
that at some particular value of s the corresponding values of {Bz,w) and {By,v) are not independent. 

It is the purpose of this paper to construct a perturbation expansion based on the smallness of e 
which allows one to put the adiabatic approximation in conte xt and t o allow for rapid oscillations on 
the small scale. These effects have been studied analytically in|ll 'ichevi lll996 h and later numerically by 
l^akholdin and Il'ichev (1991) and iBakholdi n et al... (2002.1 . However, those studies were only carried out 
for the case V — l + /i/2 where ^ is small. In the present paper, the analytic theory based on the smallness 
of £, a naturally small parameter, clearly shows the origin of all the basic features of the solutions they 
obtain and is not restricted to small fi. 

The underlying form of the governing equations suggests that a multiple-scale perturbation expansion 
is appropriate and this is carried out in Sec. 2 where explicit equations are obtained which describe the 
evolution of v, w, By, and B^ on both the small and large scale to lowest significant order in e. This 
analysis is only valid when cos 9 is not small. For 6 close to ±7r/2 singularities develop in l|4a|l and so 
an alternative set of variables and equations derived from Hla(l must be employed, as is shown in Sec. 3. 
The final section summarizes our results and the further applicability of our approach is discussed. 



2 Multiple-scale perturbation expansion 

The perturbation expansion is based on the implication from equations (|4a|l that two distinct spatial 
scales exist - a large one and a small one. We formally proceed by introducing multiple scales, s, si — 
es, S2 = e^s, . . . and treating them as independent variables so that 

d__d_ _d_ 2 

ds ds dsi ds2 ' 
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and allowing all dependent variables / to be functions of s, si, . . . and expressible in the form 

f{s) = /o(s,si, . . .) +£/i(s,si, ...) + •••, 

although in the following it is not necessary to consider the scaled variables beyond s (the small scale) 
and si ( the large scale). For a rn ore general introduct ion to this form of perturbation theory see, for 
example, iNavfeh and Mood l|l979l) and lRowland^ l)l990l) . 

Substituting the above form for the dependent variables into H4a() and Ij4b|l gives to lowest order 

^^0 _ „ dwQ _ 
as OS 

and so vq and wq can only be functions of si, and not s. At lowest order, |0c|* and lj4d|l reduce to 

^ = no{Vwo + B.o}, ^ = no{Vvo + Byo}. (6) 
Since vo and wq are independent of s, the solution of these equations is 

ByQ = hy-Vvo, Bzo = h^~Vwo, (7) 

where 

dhy dh^ 

= rio/i^, = -nohy. 

These equations have the solution 

hy = h{s i) sin 4>, — h(si) cos4>, (8) 

where 

^ = -0 (9) 
OS 

and h{si) is an as yet to be determined function of si. Substituting the above results into ^ gives, to 
this approximation, 

— = A + D{si)cos<j) + E{si)smcl) (10) 
no 

^ 2Vta.nevo-V^iwl+v^o)-h^ ^^^^ 

VhwQ {Vvo-sme)h 
^ = ^- ^2 • 

Thus A, D, and E are functions of si only and the s variation of uq is through (/) only. Integrating the 
reciprocal of (PJ after using 1)10(1 to express uq in terms of cj) yields 

A(p + D sin (t) - E cos cj) = s - s{si) (12) 

where s is a function of si. A solution of p2|l in the form = 0(s, si) is obtained in Appendix IbI 
We now proceed to next order and find using (|la|l that 

— — = — nowo + V^wq — Vh cos (j). (13) 

OS as I 

Since the right-hand side depends on s via (f) only, we replace dvi/ds by nodvi/ d4>. Remembering that 
Vq, Wq and h are independent of cj), and using (|10|l . (|13l) is readily integrated to give 

vi = ( Wo ~ '^r^i (A(j) + D sin cj) — E cos (j)) — wo4> 

+ Vh ^ jCos2(/)- ^(2(/) + sin2(/)) - Asin^j + w(si) (14) 



where 



and 
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where v{si) is an undetermined function. To ensure that vi is a bounded function of s, we must remove 
terms proportional to 0. After replacing D by its definition this results in the consistency condition 

^77^= ^y^-i-^No- (15) 



dsi V 

This is the equation for the variation of vq on the slowly varying scale, s\. Using IjlSI) to simplify (|14|l 
leaves us with 

ui = -i- ^wq + (-D sin — cos 0) + V"/i cos 2(/) — ^ sin 20 — A sin 0^ (16) 



in which v has been absorbed into vq. 
Similarly, to first order, H4b|l yields 



O 1 

^ = — ^ V-^fo + ^/isin^ + ytang (1 - tiq)- (17) 

as dsi 



After again replacing djds by 11^8/841 and integrating, to obtain a bounded Wi we require that 



with the result that 



1 „ VhE\ 



Wi — -^yV tanO — vo ^ — 1 (13 sin — £' cos ( 



- y/i ^^cos0+ ^cos20+ jsin20^ . (19) 

To lowest order, the adiabatic approximation |(SJ) is 

Bm = ~Vwo, ByQ = -VvQ. (20) 

Comparison with equations iQ shows that this approximation is equivalent to setting h ~ 0. Equations 
(I15|l and (|18|l then form a complete set which can be integrated. The details are given in Appendix IXI 
In particular, the existence of solitary pulses is proven. 

If the adiabatic approximation is not made, it is necessary to obtain an equation for the variation of 
h on the si scale. This is achieved by considering the equations for By and to next order in e. From 
(|^ and (|4d|l we may write, respectively. 



ds dsi 
dB^^dB^ 
ds dsi 



= no{Vwi + Bzi} + nihcoscj) 
= —no{Vvi+Byi}—nihsm( 



Since the variation of all coefficients with s is through 0, we replace d/ds by nod/dcj) and rewrite the 
above equations as 

9-B„i 1 dB,.o - nih , , 

^ + TT^ = Vwi+B,i + ^coscb, (21a) 

ocp no osi no 

dB,i 1 dB,o ,7 A nih 

— — H — = -Vvi - Byi sm(/). (21b) 

0(p no osi no 

We proceed by adding (|21a|l multiplied by sin (f> to (j21b|) multiplied by cos (j) and then integrating from 
to 27r. Insisting that Byi and B^i are periodic functions of (p means that 

( „ , sm ) = — ( Byi cos ) , ( cos <p ) — (Bzi sm 0) , 
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where (•) denotes the average as cj) varies from to 2tt. The combined equations then reduce to the 
foUowing equation for the variation of h: 

l \ Ah f.U^Vf) /sin0\ dw^o /cos0\ • ,\ / ,\\ ^oo^ 

V \ - — ( ) + - — ( ) + {wi sm 4)) - {vi cos 4>) > . (22) 



^?io/ dsi [dsi \ TiQ / dsi \ no 

Using ifTUIl . and ifTOIl we have 



1 \ / sin (t>\ E I cos D 



no I \ no / 2 ' \ no / 2 

• f.r /, VhE\ D , { V^h'^\ Ewo 

(wism0) = ( ytane'-wo — I (wiC0S(/)) - (1 + ^^ I 

Then inserting the above expressions and results (|15|l and 118|l into (|22|l and simphfying we find that 
the right-hand side of H22|l is zero and hen ce that /i is a constant. This is in agreement with the result 
obtained for V close to 1 in lll'ichevl l)l996j) . 

In summary, we have seen that to lowest order, the ion drift velocity components v and w only show 
large-scale variation. Fast periodic variation occurs at the next order of approximation, as given by (|16|l 
and H19|l . but given the smallness of e, these oscillations would be barely discernible. On the other hand, 
even to lowest order, the magnetic field components show rapid oscillations on top of the large-scale 
variation: 

Byo -Vvoi-si) + hsm(j){s, si), B^o = -Vwo{si) + hcos4>{s, si), 
where, as is shown in Appendix IbI sin cj) and cosi/) are periodic functions of the variable S given by H46|l . 



3 The small cos^ and 6 = ±7r/2 limits 

So far we have treated cos as finite but our treatment does not allow one to pass to the case cos 6 = 
since in this limit V becomes infinite. To consider this limit we write 9 = ±7r/2 =F y/eip so that cos 6* = 
sm{^/e'^p) ~ y/eip. (In the remainder of this section the upper and lower signs refer to the cases where 9 
is in the neighbourhood of tt/2 and —tt/2, respectively.) It is now necessary to go back to the original 
equations expressed in terms of ^ rather than s and define X = \/RiRe^. Also, to avoid singular solutions 
we need to use scaled versions of the ion drift velocities, defined by u = v/\/£ and w = wj^fe. Then 
dividing Hla() by ^/RiRl we obtain 

% = -S, + 0(£V), (23a) 



By±l-n-f-0(eV ), (23b) 



AX 
dw 

'ax 

^ = nw+'^^l^ + Oie^^), (23c) 
_ nBy 3 



AX V 



?/'-hO(eV ) (23d) 



with 

1 1 



^^=l-:^[Bl±2By + Bl)+0{e4'^). (24) 

Any solution to the above equations will be such that v ^ B and hence the real velocity v = yfev will be 
small compared to B. 

For the 9 = ±7r/2 limits, equations (|23a|l reduce to 

-Bz. = -nv, (25a) 



AX " AX 

Al 

= By±l-n, — 
AX ^ 'AX 



"^"^ By±l-n, ^^nw. (25b) 



We now define an operator L by 



AX \nAX 
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Then from H25al) and Ij25b|l it can be seen that, respectively, 

LB, = B„ (26a) 
LBy = By-n, (26b) 

in which we have re-instated By, as given in (j^J. Combining equations H26a(l gives 

ByLB, — B^LBy = B^fi. 
Integrating this over one period (or all X if boundary conditions permit) yields 

{B,n) = 0, (27) 
where (•) denotes the integral over X. Similarly, H26a|) and l|26b|l imply that 

(S,) - 0, {By) = (n) . (28) 

Relations (|27|l and H28(l are satisfied if n and By are even functions and B, is an odd function of X. 

We can demonstrate the existence of a non-trivial solution of (|25a|l by taking B, = 0. In this case 
(|^ becomes 

- = a{l~ PBl) 
n " 

where 

with the result that Ij26bp can be re-expressed as 

"dx((^-^^^)df)^^^-;^(T^- 

Letting p = dBy/dX, this can be written as the following first-order differential equation for p^: 



dp2 / 4:f3By \ 2 _ 2 / B. 



dBy \1-[3BIJ^ a\l~[3Bl a{l - PB^)^ 

32\2 



which has an integrating factor of (1 — (3By) . Hence the solution is given by 

= n ^^Rm,2 (30) 

where 



B,, 0B,i B,, 

Q = ^--^-^ + Qo (31) 

and Qo is an integration constant. For a solitary pulse solution, the appropriate boundary conditions are 
By ^ ±1 and Q ^ as |X| — > oo. Using these allows us to determine Qq. We then rewrite H31() as the 
following expansion in By: 

Q = Bl[^Tf3By-^Bl]. (32) 

A necessary condition for the existence of solitary pulses is therefore that f3 < 1/3. In addition, from (|30|l it 
can be seen that p is singular at By — 1/ (3. Hence, for a solitary pulse to exist, at least one of the two non- 
trivial zeros of the expression for Q given in 1)321) must lie within the range — 1/v^T ^ < By < l/vT^T 1- 
The zero which is larger in magnitude never satisfies this. The remaining zero at ±v/2//3- 2 =F 2 satisfies 
the condition if /3 > 1/9. Using (|29)l we can now write the sufficient condition for the existence of a 
solitary pulse solution as 1 < U^ < 4. 

As a check on our calculation we look at the case when V = 1 + /i/2 for small positive ^. Then (|3()|l 
reduces to 

to lowest order in /i. The above equation has the solution 

By = i^sech^ \^/T^{X - X^) + 0(^/2) 
where Xq is an arbitrary constant. This is in agreement with the result given in lBakholdin and irichevl 
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4 Conclusions 



We have studied a set of magnetohydrodynamic equations for planar magnetoacoustic waves of permanent 
form propagating in a two-component cold plasma and, by taking advantage of the smallness of the 
ratio of the electron to ion masses, have obtained a reduced set of equations which describe the large- 
scale variation of the magnetoacoustic wave solution of the full equations. Superimposed on the large- 
scale variation, multiple-scale perturbation analysis indicates that there is a rapid oscillation which is of 
constant amplitud e in the c ase of the lowest-order magnetic field components. These results are consistent 
with the study of lll'lchevl l(l99ft l which was restricted to a narrow range of velocities. In addition, the 
approach expounded in this paper puts the adiabatic approximation into its true context. 

In this work we have obtained various conditions for the existence of solitary pulses. Whether these 
solutions correspond to phenomena that could occur in nature depends on whether they are stable. Linear 
stability analysis of the solutions shown to exist in this paper is a challenging problem. However, the 
numer ical solution of the full (time-dependent) system of equations obtained bv lBakholdin and 11' ichevi 
ljl998l) . show that for a range of initial conditions the solution relaxes to the type of solution shown to 
exist here. This suggests that our solutions are stable, at least to perturbations applied in the direction 
of propagation. 

Although the equations studied here arise from a magnetohydrodynamics problem, the method is 
applicable to a more general set of nonlinear equations where two distinct scales are a basic feature. 
An advantage of the present study is that the equations obtained on the large scale can be investigated 
analytically and describe real physical processes. Although we have only looked at cold plasmas, an 
exactly analogous procedure can be applied to the case of warm plas mas, at the expens e of som e additional 
algebraic complexity. The relevant governing equations are given in lBakholdin et ahl l)2002(l . 
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A Adiabatic approximation 



When h = 0, the coupled equations for vq and wq, namely 115(1 and ((18|l . reduce to 




(33) 



and 




(34) 



Multiplying and by 2vq and 2wo, respectively, and adding gives 




(35) 



Differentiating Hll|) with respect to si and using the above, one obtains 




(36) 



(37) 



and 



(38) 
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Figure 1: Phase plane for H4U|) when conditions are such that both compressive and rarefactive sohtary 
pulses occur. 



where Ci and C2 are integration constants. Finally, after combining H3t)|) . (|37|) and (|38|l one obtains 

[^^^-t.^-^- (39) 

^ ^ m— 

in which 

7o = (Ci-C|)k^ 71 = 2(C2 - ytan6l)K, 72 = tan 6I - CaF^)/* - 1, 

73 = 74 = -1^^ 

This is of the form of the energy equation of a particle with position A in a Sagdeev potential which is 
minus the right-hand side of H39|) . In general, nonlinear waves exist and in particular solitary pulses. The 
latter can occur when the boundary conditions are such that u — > 0, w ^ 0, and n — > 1 as si — > ±00. 
Using the result that in this limit A ^ 1, the boundary conditions allow us to use H37|l and H38|l to 
determine the integration constants in this case: 



Ci=VK-HaTi0, C2 = [l- ^V^)k~^. 



Then reduces to 



A—] = ^ {a + bA + cA') (40) 

dsi / 



where 

The requirement that 14U|) gives rise to a solitary pulse is that a + 6 + c > 0. Using the above renders 
this condition as (V^^ — 1)(1 — V^) > which on rearranging yields 

cos^e* < < 1. 

As illustrated in Fig. ^ since c < 0, if the above condition is satisfied, compressive solitary pulses will 
always occur. However, rarefactive pulses are only possible if the smaller root oi a + bA + cA? = is 
above zero. This will occur if a < and V"^ < 2. These requirements are equivalent to the condition 

< 2(1- I sin 61 1). 

This implies that if the values of V and 9 are such that compressive pulses exist, then rarefactive pulses 
will also occur if 16*1 < 7r/6. 

It is possible to integrate H4U|) to obtain the spatial variation of the solitary pulses implicitly. An 
approximate explicit solution can be obtained when is just above 1. Introducing U = A — 1, (|4U|I 
becomes 

'duV V^U^{U+-U){U -U-) 



dsij (1 + C/) 



2 



(41) 



where 



t/± = ?^^^(±tan^-y^^ 



8 



If - 1 « I tan 6* I, then U± ~ ±v where 



2VT723Ttan6l 



y2 

For soUtary pulse solutions, \U\ < \U±\^ and so if tan 6* is of order unity, C/«l. Hence (|41|l reduces to 



dsi 

at lowest order and one obtains 

A~l ±vsechvsi. (42) 

Using H36|l and (|38|l we can then obtain the corresponding expressions for wq and vq: 

Wo — T4V^tan6'sechzysi tanhz^si, wq — 2V^tan0sech^ I'Si. (43) 

In this adiabatic approximation, the lowest order components of the magnetic field are just multiples of 
these quantities, as given by 12U|) . 

We can now also use (|42ll and H43() to obtain the solution when — 1. This corresponds to the limit 
— > in which case A ^ 1, — > 0, and wq ^ 2Vi&n9. As a check on our calculation, we note that 
these results are consistent with the definition of A as given by when /i = 0. 

B Explicit expression for 0(s, si) 

The variation of (j) with s is given implicitly by (|12|l . An explicit expression can be obtained by writing 
the equation in the form 

S + ^] = (l) + iP-<7 sin(0 + V) (44) 



where S = {s — s)/A, = arg(— J + i£') and a = V + E'^/A. The following explicit solution to 
was first obtained by Ijacksoni l)l960|) (although for a more transparent exposition see p. 154 of 



El 

llnfc: 



llnfeld and Rowlands (,2000,) ): 

(/) = 5 + 2 y '^"^'''^^ sinm(5 + V'). (45) 

m—l 

It can be seen that (j) has a directed component, S, on which a periodic variation is superimposed. Since 
A and hence the period vary on the si timescale, to the order in e to which (|45|1 applies, it is more 
appropriate to re-define S by 

Such a definition avoids secular terms at higher order in the e expansion. 

References 

Bakholdin, I. and Il'ichev, A. 1998 J. Plasma Phys. 60, 569. 
Bakholdin, I., Il'ichev, A. and Zharkov, A. 2002 J. Plasma Phys. 67, 1. 
Haken, H. 1983 Synergetics. 3rd edn. Springer- Verlag, Berlin. 
Il'ichev, A. 1996 J. Plasma Phys. 55, 181. 

Infeld, E. and Rowlands, G. 2000 Nonlinear Waves, Solitons and Chaos. 2nd edn. Cambridge: Cambridge 
University Press. 

Jackson, E. A. 1960 Phys. Fluids 3, 831. 

Kakutani, T., Kawahara, T. and Taniuti, T. 1967 J. Phys. Soc. Japan 23, 1138. 

Kakutani, T., Ono, H., Taniuti, T. and Wei, C.-C. 1968 J. Phys. Soc. Japan 24, 1159. 

Nayfeh, A. H. and Mook, D. T. 1979 Nonlinear Oscillations. Wiley, New York. 

Rowlands, G. 1990 Non-Linear Phenomena in Science and Engineering. Ellis Horwood, London. 



9 



